Two-Qubits Entanglement Induced by a Faster Data Bus 
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We propose a theoretical protocol to create the entanglement of two qubits via the Born- 
Oppenheimer (BO) approximation. In our scheme, each qubit is coupled to a faster data bus whose 
frequency is much larger than the energy spacing of the qubits and thus the BO approximation is 
valid. Then the adiabatic separation of qubits from the data bus can induce an effective potential 
to couple the two qubits, which can be utilized to create a quantum logic gate. We also discuss 
the quantum decoherence caused by the adiabatic entanglement between the two qubits and the 
external field. 
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I. INTRODUCTION 

In quantum information theory, it is important to cre- 
ate a "maximum entanglement" between two qubits. It 
has been shown by Barenco et. al. [1] that, together 
with single-bit operations, the entanglement based two 
qubit gates except the classical SWAP gate (including, 
the CNOT [2] and controlled phase gates) forms a univer- 
sal set of logic gates of quantum computation. In other 
words, any unitary transformation of single or many- 
qubit system can be decomposed into the quantum net- 
works of single-bit gates and non-trivial two-qubit gates 
. In the last few years, many schemes to create two-qubit 
interaction have been suggested and experimentally im- 
plemented [3-8]. In this paper, we will propose a new 
scenario to induce a kind of two qubit interaction cre- 
ating efficient quantum entanglement. In our proposal, 
each qubit interacts with a common quantum field as a 
data bus. If the frequency of the quantum field is much 
larger than the energy spacing of the qubit, the degree of 
freedom of the quantum held and the variables of the two 
qubit system can be separated adiabatically. In this sit- 
uation, an effective two-qubit interaction can be induced 
by means of the Born-Oppenheimer (BO) approximation. 
To explain the main idea in our proposal in details, we 
will hrst briefly review the relevant quantum adiabatic 
theorem and the BO approximation as follows. 

Adiabatic approximation [9] is widely used in many 
fields of physics including classical physics and quan- 
tum physics. In the latter, one has a rich applications 
of quantum adiabatic approximation (QAA). We con- 
sider a time-dependent quantum system. Its Hamilto- 
nian H (t) = H [B4 (t)] is controlled by some adiabati- 
cally changing parameters Ri (t)s. If the system is ini- 
tially prepared in the n-th instantaneous eigen-state of 
the Hamiltonian H (0) at t = 0, then during the evo- 
lution the quantum system will still persist in the n-th 
eigen-state of the instantaneous Hamiltonian H (t) at an- 
other instance. The evolution of the system under the 
control of Ri (t)s can be evaluated by QAA, but the back 
action of the quantum system on the controlling param- 



eters Ri (t) is neglected. Of course, the QAA can only 
work well when Ri (t) behave classically. 

When the Hamiltonian of the quantum system depends 
on the variables of another slowly varying quantum sys- 
tem, (e.g., the motion of an electron near nuclear), one 
needs a more exact approximation, in which the inter- 
actions between the fast system and the slow one, espe- 
cially for the "back action" of the fast one on the slow 
one, are considered. To this end, in 1930 [10], Born and 
Oppenheimer suggested a high order approximation (we 
now call it Born-Oppcnhcimcr (BO) approximation) to 
tackle this problem. In the view of the BO approxima- 
tion, to solve the dynamics of a composite system with 
fast and slowly changing parts, we can hrst solve the 
Schrodinger equation of the fast part for each fixed slow 
variable q. Then the obtained g-dependent eigen-value 
V n (q) provides the slow part with an effective potential 
when the motion of the slow part can not excite the tran- 
sitions among those fast states, i.e., the adiabatic condi- 
tion holds. In this sense, the time dependent wave func- 
tion of the composite system can be written as 



V(q,x;t) = ip n (q,x)<j>(q,t). 



(1) 



Here, x is the variable of the fast part, ip n (q,x) the 
q-dependent eigen-state of the fast part corresponding 
to eigenvalue V n (q); and the slow part <p(q,t) evolves 
governed by the effective Hamiltonian Hf +V n (q), i.e. 

cj)(q,t) = e - 4 [^ +y " (,5) ]V((7,0), where H s is the free 
Hamiltonian of the slow part. 

In this paper we make a crucial observation for the cre- 
ation of the nontrivial two qubits interaction of logical 
gate in quantum information as a result from the appli- 
cation of the BO approximaiton to a three body system. 
According to above discussion, it is very interesting that 
the slow part can be divided into two subsystems with 
variables q\ and q 2 without interaction between them. In 
this situation, the Hamiltonian of the composite system 
can be generally written as 



H = H 1 {qx) + H 2 {q 2 )+W l {q u x) 



W 2 (q 2 ,x)+H f (x) 
(2) 
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where Hi (qi), H% (q 2 ) and Hf (x) are the free Hamilto- 
nians of the two slow subsystems ( we can note them as 
(71,(72 ) and the fast subsystem. Wi{q\,x) (W 2 ((72, %)) 
are the interaction Hamiltonians between qi ((72) and x. 
As we have emphasized, Eq. (2) implies that there is 
no direct interaction between q\ and qi. Under adia- 
batic condition, the eigenvalue V n {q\,q2) of the Hamil- 
tonian W\ (q\,x) + W 2 {q 2l x) + Hf (x) of fast part in- 
duces an effective potential on the slow variables qi and 
(72. In usual case, V n (qi,q 2 ) can not be decomposed into 
the sum of two independent potentials of qi and (72 i.e. 
V n (qi,q 2 ) ^ Vi n (qi) + V 2n (q 2 ). Therefore, an effective 
interaction between q\ and q 2 is induced via the BO ap- 
proximation and then one can regard the fast variable x 
as a data bus which can be removed. Indeed, in quantum 
information process, the inter-qubit interaction is very 
crucial to create the nontrivial interaction between two 
qubits q\ and q 2 if we can regard the above mentioned 
two subsystems of q\ and q 2 as two qubits. 

We should point out that, in many previous scenar- 
ios of quantum logical gate, an external field was usually 
used as a data bus where the quantum information car- 
ried by a qubit is stored and then transferred to another 
qubit. Most of those proposals rely on the rotation wave 
approximation (RWA) . The effective interaction between 
the two qubit system can be obtained through the adia- 
batic elimination of the external field [8]. 

However, to adiabatically eliminate the degree of free- 
dom of the external field, the coupling strength g be- 
tween the qubit and the external field ought to be weak 
enough so that the ratio £ is very small. Here, A is the 
detuning i.e. A = \uj a — lo\ where ui a is the transition 
frequency of the qubit and lo is the frequency of the ex- 
ternal field. On the other hand, since the RWA requires 
A < Lu a , the above large detuning condition implies that 
g << oj a - Then the effective coupling strength between 
the two qubit system ^ is much smaller than uj a . This 
weak coupling may reduce the "quality factor" Q of the 
logic gate. Here Q can be understood as the ratio be- 
tween the atomic life time and the logic operation time 

77 ^75^ since the Rabi transition frequency can be 
defined as an upper limit of the operation frequency. 

In our proposal based on the BO approximaiton , the 
frequency of the external field can be significantly larger 
than the energy spacing of the qubits. Then g may 
be larger than uj a or, at least, have the value in the 

same order as uj a . The ratio between effective interaction 
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strength and u> a can be much stronger than that in the 
previous scenarios mentioned above. In the next section, 
there is a brief introduction of our quantum logic gate 
model. The precise analysis of our model based on the 
perturbation theory is given in section III. In section IV, 
the adiabatic decoherence effect is analyzed detailedly in 
association with the idea of the adiabatic entanglements 
presented recently by Sun et. al. [11]. We also obtain the 



condition under which the decoherence effect can be ne- 
glected and the two-qubit interaction model works well. 
There are some conclusions and discussions about the 
physical realization of our protocol in the last section. 

II. EFFECTIVE INTERACTION BETWEEN TWO 
QUBITS 

We generally consider a system consisting of two qubits 
with the same energy spacing u a and an external field 
which is described as a dynamic data bus. We can treat 
the external field as a harmonic oscillator with frequency 
w (see Fig. 1). 



Qubitl Data Bus Qubit2 




FIG. 1. The interaction between two qubits and a external 
field. When uj >> u a , the external works as a faster data bus 

The Hamiltonian of the composite system reads 

H = uj a (of + of) + g (fit + a) (of + of) + u tfa 

(3) 

where of = |e). (e| - \g) l (g\, of = |e). (g\ + \g). (e| , |e), 
and \g) i are the excited and ground states of the i-th 
qubit; a? and a the creation and annihilation operator of 
the external field, and g the coupling strength between 
the qubits and this field. This model system can be real- 
ized in many practical physical situations such as the two 
trapped ion with the spatial oscillation mode, two atoms 
in a single mode cavity and the two Josephson charge 
qubits coupling to a large Josephson or nonamechanical 
resonator [12] . We assume u> is much larger than uj a . In 
this case, the RWA could not be applied since the con- 
diton for the validity of RWA is \u> + uj a \ >> \u — u) a \. 
Thus, the usual adiabatic elimination technique does not 
work well for this case. We note that the above Hamil- 
tonian depicts a typical spin-boson (continuous variable) 
system. 

In the following parts of this paper, with the gener- 
alized BO approximation [11], and perturbation theory, 
we will prove that under the conditions lu » cu a and 
lu » g, the evolution of the two qubit system is approx- 
imately governed by the effective Hamiltonian 

h = Ua(*P+*W)-2£*W*j?> (4) 
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and the decoherence effect caused by adiabatic entangle- 
ment [13] between the two qubits system and the data 
bus can be neglected. 

Under the adiabatic condition that the external field 
does not excite the transitions among the internal states 
of qubits, the BO approximaiton determines a formally 
-factorized eigen-state of H for the total composite sys- 
tem. 

(Ai,A 2 I*} = <p na (Ai,A 2 ) |n[Ai, A 2 ]) (5) 

where |Ai,A 2 ) is the common eigen-state |Ai,A 2 ) 
of |o-i 1 ' ) , a x 2 ^ | with eigen-values {Ai,A 2 =±l} and 
|n[Ai, A 2 ]) the eigen-state of the effective Hamiltonian 

Hf[X u A 2 ] = Lua^a + g (a f + a) (Ai + A 2 ) (6) 

of the external field for the given state | Ai , A 2 ) of the two 
qubit system with corresponding eigen-value 

V n (XuX 2 ) ^ ncu - ^ (X 1 + X 2 ) 2 (7) 
g 2 

= nuj — 2 : — AiA 2 . 

Here, we have neglected a constant term. For different 
n, V n (Ai,A 2 ) contribute different effective potentials for 
the slow part so that the wave function <j> na of the two 
qubit system is determined by 

H eff(n) \4> na ) = E na \4> na ) (8) 

where the effective Hamiltonian is 

H en {n)^H eft in-a^\a^) (9) 

= . a (^)+a( 2 ))+K(4 1) ,-i 2) ) 

= h + nui 

where h defined in Eq. (4) is the n-indepedent part of 
the effective Hamiltonian. 

So far, we have formally obtained a X — AT-type of 
interbit coupling based on the BO approximation, which 
can be conveniently used to formulate an efficient scheme 
for quantum-computing. Recently, a similar interaction 
generating the entanglement of two Josephson junction 
charge qubits was proposed by Averin [14]. In ref. [14], 
two Josephson junction charge qubits are assumed to be 
coupled with another junction which works as a faster 
data bus. With the BO approximation, the energy of 
the lowest band of the latter junction can be considered 
as the effective interaction between the two Josephson 
junction charge qubits. 

In the next section, we give a detailed derivation of 
the effective Hamiltonian (9) based on the perturbation 
theory. The decoherence of the quantum state of the two 
qubit system induced by the adiabatic entanglement with 
the external field is analyzed in section IV. 



III. GENERALIZED BORN-OPPENHEIMER 
APPROXIMATION FOR DISCRETE SYSTEM 



In the above section, we obtain the effective Hamil- 
tonian (4) for two qubit system through a intutitional 
analysis from the BO approximation directly. However, 
the BO approximation is usually used to seperate two 
continuous quantum systems (e.g. a melocular and a 
electron) or a continuous system and a discrete one (e.g. 
the spatial motion and spin of a neutral particle [13]), 
rather than the two descrete subsystems. In this section, 
however, we can generalize the BO approximaiton to deal 
with the seperation of two discrete quantum systems, i.e. 
the two qubit system and the external field. We now can 
give a rigorous proof to the intuituional argument in the 
above section based on the perturbation theory. In this 
way, we obtain the effective Hamiltonian (9) naturally. 

A. Matrix Representation of Motion Equation 

In section II, we define |n(Ai,A 2 )) as the eigen-state 
of the Hamiltonian _ff/[Ai,A 2 ] of the external field and 

|Ai,A 2 ) the eigen-states of jai^, j. With the com- 
pleteness relation [11] 

E E \n (Ai, A 2 )) (n (X U X 2 ) | ® |Ai, A 2 ) (A 1; A 2 | = 1, 

Ai,A 2 n 

(10) 

the quantum state of the composite system can be ex- 
panded as 

|*>= E C Xl , X2 (n)\n[X 1 ,X 2 }) (11) 

{n,X u X 2 } 

where 

\n [Ai, A 2 ]) = |n (Ai, A 2 )) |Ai, A 2 ) . 

Substituting Eq. (11) into the Schrodingcr equation 
H\*i?) — E\ty), we obtain the effective equations of the 
coefficients C\ lt \ 2 (n): 

E H K'il (») C K-M (») + E F K,K (») C KK (*) ( 12 ) 
+ E E °x{il ^ m ) C K,K (™) - EC ^M (n) 

X[ ,X' 2 m^n 

where we have defined some complicated notations: 

HxiH (») = <^i, A 2 | H eff (n) \X[,X' 2 ) (13) 

Fx>£ H = (Ai, A 2 | w« (*P + of ] ) | Ai, A 2 ') (14) 
x[l-(n(A 1; A 2 ) |n(Ai,Ai))] 
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O x x }> x x ? (n,m) - (Ai,A 2 |w„ (<r« +^ 2) ) |Ai,A 2 ) (15) 
x (n(Ai,A 2 ) |m(AU)). 

Obviously, only when F x , ' x ? (ri) and Oyp(n, m) 
are negligible, can we obtain the effective Hamilto- 
nian H e ff(n). With a straightforward calculation (see 

A 1 A 

the Appendix A), (n(Ai,A 2 ) \m(X\,X 2 )), F x , ' x , 2 (n) and 

O x }\? (n, m) can be expressed as power series of the pa- 
rameter — : 



(m(A 1 ,A 2 ) \n(X[,X' 2 )) 



(16) 



exp 



1 




A 2 


2" 


~2 




A 2 





E 



J!(/ + n-m)! 



2i + (n-m) 



\/n (n — I) ... (m — I + l)\/m (m — 1) .. (m — I + 1), 



— \m— E v \, o 



n.m±l 



(17) 



(18) 



Ai + 



where we have assumed n — m > and Ew'w 
A 2 - A' x - A' 2 . According to Eq. (17), F^ (n) can be 

neglected if n and — are small enough. Considering the 
condition u> a << lu mentioned in the above section, in 
the following discussion, we assume n < 1 , g << u> and 
o> a << w. For simplicity, we also assume g~ui a and 
is a small parameter. 

In order to solve Eq. (12) with the standard perturba- 
tion theory and thus give a proof for the generalized BO 
approximation for the discrete case, we rewrite Eq. (12) 
in matrix-value form as 



g ~Wg 



EC = EC 



(19) 



where the total Hamiltonian E and the eigen-vector C 
are defined as 



E 



H(0) + F(p), 0(0,1) 

6(1,0), S(i) + F(i) 



,c 



Ci 

c 2 



Here, the vector C contains two sub-vectors 



d = ( Cu W , Ci_i (i) , C-u (*) , C-i-i (i) ) , i = 0, 1. 

Each of them is a 4-dimension vector; F (i) and 

6 (i, i) (i, j = 0, 1) arc 4x4 matrices whose elements are 

Making use of Eq. (17) and Eq. (18), we express the 
matrix H as a power series of -: 



where 



H (0) 



0, 
0, / 



o, 



and 



0(2) 



~(l)-(2) 



ct (1) «t (2) 



9 ' 



(20) 



(21) 



, (22) 



(23) 



(24) 



Here, / is the unit 4x4 matrix and E a 4 x 4 matrix whose 
elements are E^J' A ?. Only the terms of zcroth, first and 

second order in — are considered in Eq. (20). The terms 
of higher order are neglected. According to Eq. (17), 
F (i) (i = 0, 1) docs not contain the terms of zeroth, first 
and second order. Therefore, F (i) are neglected in Eq. 
(20). 



B. Perturbation Solution of Motion Equation 

To apply the perturbation theory, we split E into two 
parts: 



H = E + W 



(25) 



where 



Ho = + £ (26) 

is treated as the unperturbated Hamiltonian and 

W=(£f (V^+O^) (27) 

as the perturbation term. 

It is pointed out that we can also choose E^ as the 
unperturbated Hamiltonian, ^ijW the first order per- 
turbation term and (^) 2 (V^ + O^) the second order 
one. It is easy to prove that, if B-H^ is regarded as 
the first order perturbation Hamiltonian, the eigenenergy 
and eigenstate to the first order given by the perturbation 
theory is just the same as the ones we obtain by diago- 
nalizing the matrix E^ + jjjE^ exactly. Therefore, our 
final result is not relevant to the choice of unperturbated 
Hamiltonian and perturbation terms. 



4 



The eigen-energy and eigen-state of unperturbated 
Hamiltonian Hq in Eq. (26) are easy to obtain. The 
eigen-energy can be expressed as: 



E {0) -e- 



(28) 



where n = 0, 1 and i = 0, +1, — 1. Here, we have defined 

£o=0, e ±1 =±2iu a . (29) 

The energy level structure of the compositive system 
is illustrated in Fig. 2. Since to a is much smaller 
than u>, it follows from Eq. (28) and Eq. (29) that 



E 



(0) 



E, 



(0) 

oj 



uj. Notine 



E 



(0) 



we have 



E 



(0) 



(0) 



■n i 
,(0) 



(0) 



(0) 



Therefore, the en- 



<< E u -^ 0j 

ergy spectrum of Hamiltonian possesses a qusi-band 
structure: it consists of two qusi-bands corresponding to 
n = and 1. 



B = 



FIG. 2. The energy levels of H^K Since ui >> oj a , we have 
— ^0j' >> £j ~ £j- Therefore, there are two subbands 
with respect to n = and n — 1. 

The eigen-states corresponding to iJ n ± 



4°2i) =(c ± , Of. 



4%) = ( 0, C± ) T (30) 
can be written in terms of the unit vectors 

C+ = (l, o, o, o) T , 

C7_ = (0, 0, 0, 1) T . 

It is obvious from Eq. (21) and Eq. (22) that the eigen- 
energies E^q are two-fold degenerate, which corresponds 
to the two degenerate eigen- vectors 

\\E 00;k ) =(C k , 0) T , 
ll^iojfc) =(0, C k ) T . (31) 
where k = 1 , 2 and 



d = ( 0, 1, o, o ) T 
C7 2 = (0, 0, 1, 0) T . 

By making use of the perturbation theory, the eigen- 
energies and eigen-states to the first order can be written 
as 



E - E {0) 



E 



(1) 



and 



ll^nil) 

\E n Q- : k) 



E 



E 



(0) 
nil 

(0) 



E 



E 



(i) 

n±l 

(1) 
i0:k 



According to the perturbation theory for degenerate lev- 
els, the zero-th order eigen-states 
binations of \\E n Q-i) and | E n , 

(£„ ;l||O( 2 )||£„ 0;2 ) = 0, \E t 

diagonalizing the 2x2 matrices 



■0;2/ 
(0) 



^no-k/ are l mear com- 
For any k and k', 
can be obtained by 



(E„ -i\\ VW \\E n0 -i) , (E n0 ;l\\ \\E n0 . 2 
(E n o ;2 \\V^\\EnO;i), (E n0 -M\V^\\E 



,i():2) 



(32) 



The first order corrections to the eigen-energies E n { 
and eigen-states can be calculated directly. Since we have 
(J5 n ±i|| \\E n ±i) = for every given n, it is obviously 
that the contributions from to the correction of the 
eigen-energies are zero. On another hand, the first order 

eigen- vector E^) \ can be expressed as: 



E 



(i) 



E 



(0) 



E 



(o)\ 



fe=l,2 

E 



El 



(0) 
n0;k 



E 



(1) 



(33) 



y( 2 



E 



(0) 



/^(O) 


0(2) 


^) 




(-l) n+1 u 

/ P (0) 
X^ln— l|0;fc 


+ e 4 - 

0(2) 





,(1) 

i0;fc 



(i) 

\n-l\j 



M) 



n+1 



W + Si 



E, 



(i) 

l|0;fe 



fe=i,: 

The another eigenstate End-k) nave the similar expres- 
sion and thus for simplicity, we need not give its explicit 
expression here. 

Substituting Eq. (29-32), into Eq. (33), we can cal- 
culate the contribution from and to the first 
order correction of eigen-state. It is obvious that 



E 



(0) 



E 



(0)^ 



5 



(5 
(5)" 



E 



(0) 
i0;fe 



E 



(0)\ 
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(34) 



and 



'2-X 



E 



(o) 

\n-l\j 



0(2) 



(0) 



(-1) 



n+1 



(0) 

|n-l|0,fc 



0(2) 



(0) 



(-1) 



n+l 



(35) 



where we have used — ~1. Therefore, it can be seen from 



Eq. (33) that the contribution from to 



E 



(i) 



and 



End k) are mucrl smaller than the ones from V"' 2 ' 



In summary, we have shown that one can reasonably 
neglect the first order corrections from to both the 
eigen-values and the eigen-states. 



C. The Generalized BO Approximation From 
Perturbation Theory 

The above result can be obtained with intutitve anal- 
ysis. It can be seen from Eq. (24) and Eq. (23) that the 
matrix leads to the transitions in the same energy 
band and O' 2 -* leads to the interband transitions (see Fig. 
3). 



V 1 



n = 1 



O 



FIG. 3. The transitions induced by W. We have shown 
that, (•£■) <5' 2 ' only induces the transitions between subbands 

and (f;) 2 only induces the transitions in the same sub- 
band. 

As we have mentioned, the energy gap between the two 
bands has the same order of u and is much larger than 
the energy difference between the levels in the same band. 



The latter have the same order of ui a . Then the probabil- 
ities of the interband transitions are much smaller than 
the ones of the transitions in the same energy band. 

With the above discussions, the Hamiltonian matrix 
H can be written as 



(36) 



Which is just the effective Hamiltonian for the BO ap- 
proximation. In fact, substituting Eqs. (23-21) into Eq. 
(36), we can write the eigen equation (19) as 



C = EC. 



This equation equals to 

]T H*$? (n) C K>K (n) = EC XlM (n) 



(37) 



A. ,A 



which is approximated form of Eq. (12). It can be seen 
from Eq. (13) that Hy'y (n) are the matrix elements 

of the effective Hamiltonian H e ff(n) in Eq. (9). It is 
pointed out that, Eq. (37) is equivalent with Eq. (5) and 
Eq. (8) in the above section. 

To solve the time-dependent Schrddinger equation 
which describe the time evolution of the composite sys- 
tem, we can calculate the approximate eigen-energies E 
of the Hamiltonian of the composite system with Eq. (9): 



E n 



nuj 



(38) 



where is the z-th eigen-energy of the Hamiltonian h 
defined in Eq. (4). Eq. (13) is the time-independent 
Schrodinger equation under the BO approximation. The 
time-dependent Schrodinger equation can be written as: 



ijC XuX2 (n)=£<&(«)tfAi,Ai (r 



(39) 



A. ,X 



Making use of the above BO approximation for discrete 
system, we can calculate the evolution of the composite 
system with Eq. (39). If the system is initially prepared 
in the state 

|*(0)>= Cx 1 ,x 2 (n,0)\n[X 1 ,X 2 }), 

n,(Ai,A 2 ) 

then with Eq. (38), Eq. (9) and Eq. (4), the quantum 
state at any instance can be expressed as: 



|*(t)) = J2 CAi,A a (M). 
n,(Ai ,A 2 ) 

n,(Ai,A 2 ) 



|n[Ai,A 2 ]) 



|A l! A 2 )®|n(Ai,A 2 )) (40) 



6 



where 



where 



Cai,a 2 (M) 

= E (Ai,A 2 |e^ t |A' 1 ,A^)C Vl ^(n,0). (41) 

Apparently, the evolution of the two qubit system is gov- 
erned by the Hamiltonian h in Eq. (4) which contains 
the two qubit interaction — 2^&^P &^ ■ However, it can 
be seen from Eq. (40) that because of the dependence 
of \n (Ai, A 2 )) on {Ai, A 2 }, the instantaneous state |\E f (t)) 
of the composite system is usually an entangle state of 
the two qubit system and the external field. Thus, if one 
trace over the variable of the data bus, then the the state 
of the two qubit system is not a pure sate but a mixed 
one. This adiabatic decoherence effect will be analyzed 
in next section. 

IV. DECOHERENCE EFFECT FROM THE 
ADIABATIC ENTANGLEMENT 

In the above section, we have obtained an approximate 
Hamiltonian H e ff (n) = h + nuj for the composite sys- 
tem with two qubits interacting with an external field. As 
we have mentioned, there is quantum entanglement be- 
tween the states of the external field and those of the two 
qubit system because of the dependence of |n (Ai, A 2 )) on 
{Ai,A 2 }. This is just the adiabatic quantum entangle- 
ment which has been studied in detail for the develop- 
ment of BO approximation [11]. It is easy to see that the 
adiabatic entanglement may cause the adiabatic decoher- 
ence of the quantum state of the two qubit system. Only 
when this decoherence effect can be neglected, can our 
effective Hamiltonian h be used to create a two bit gate. 
In this section, the adiabatic decoherence effect in our 
problem will be analyzed in detail. It will be shown that, 
in some practical cases, this kind of decoherence effect is 
negligible and then the faster quantum entanglement can 
be created via the BO approximation. 

We assume the two qubit system and the external field 
are initially prepared in a factorizable state: 

|*(0)) = |V>®|$ 

where is the initial state of the two qubits and \<j>) 
the state of the external field. |^ (0)) can be expanded 
with respect to the basis \n [Ai, A 2 ]) as 

|*(0))= E Ck.Ajn.OjfcHAi.Aa]) 

n,(Ai,A 2 ) 

or equivalently in the |Ai, A 2 ) picutre 

(Ai,A 2 |*(0)> =^C Al;A2 (71,0)0 1^^2]) 

n 



Cx lM (n,0) - (Ai,A 2 |V) • (n(Ai,A 2 ) \<j>) . (42) 

is the projection of inital state onto the completeness ba- 
sis |n[Ai,A 2 ]). With the BO approximation generalized 
in the above section, the evolution of the composite sys- 
tem is depicted by Eq. (40). 

Substituting Eq. (42) into Eq. (40), we have the wave 
function in the (Ai,A 2 ) picture at time t. 

(Ai,A 2 |*(i)) = 

E (A 1 ,A 2 |e-^|A' 1 ,A 2 )(A' 1 ,A 2 |V)®|^(i)) (43) 

A'i ,X' 2 

where 

\4>{t)) = |0(t,Ai,A 2 ,Ai,A£)) 

= E e'"^{n{\^\' 2 )\4>)\n{\\\ 2 )) 

n 

depends on both coordinate sets (Ai, A 2 ) and (A' 1; A 2 ). It 
should be noted that the factor (n(A' 1 , A 2 ) \<f>) \n{\\, A 2 )) 
has different pairs of indices (Ai,A 2 ) and (A' 1; A 2 ). If 
\4>{t)) = \<j> (t, Ai, A 2 , A' l5 A 2 )) were independent of Ai,A 2 
and A' 1; A 2 , then we would have a factorized evolution de- 
scribed by |* (t)) = e~ iflt \ip) ® \4>(t)). In this case, the 
two qubit system were just in the pure state controlled 
by the effective Hamiltonian h and the decoherence phe- 
nomenon does not occur. 

In the usual case, \<j>(t)) = \(f> (t, Ai, A 2 , A' 1; A 2 )) does 
depend on Ai, A 2 , \[ and A' 2 . Then the decoherence effect 
must result from the dependence |<^>(i)) to different sets 
of (Ai, A 2 , A' 1; A 2 ). In fact, we can calculate the elements 
of the reduced density matrix of the two qubit system 

x<0(t;Ai,Ai;/ii,/4) \<f> (*, Ai, A 2 ; m, n 2 )) (44) 
where the propagation functional 

^Ai,A 3i #*i,^ [|^>.*] = (Ai,A 2 |e-^*|Mi,M2> (Hl,H2\1>). 

(45) 

depend on the initial state \ip) of the external field. On 
the other hand, without decoherence, the reduced density 
matrix of the qubits is p' = e ~ tht \ip) e lht with matrix 
elements 

Px lM ,x> 2 (t) (46) 

= E E ^A 1 ,A 2W ,, l2 [|V),i]^ 1 ,A-Kx[IV')^] 
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Comparing Eq. (46) with Eq. (44), we measure the 
extent of decoherence by the so called decoherence factor 
[15]. 

F=( ( f>(t;X , 1 ,X , 2 ; ( i' 1 ,( Jl , 2 )\ ( l>(t,X u X 2 ; ( i 1 ,^)) (47) 

As pointed in ref. [15], if F = 1, there is no decoherence. 
It should be pointed out that, not only the norm of F but 
also its phase are crucially important to the description 
of quantum decoherence process. For example, if F = 
exp and 6 is a random phase with respect to the 
parameters Ai, A2, A' 1; A 2 ; m'd M2> i M 1 ; M 2 ; even if 1-^1 = 1, 
p'xt x 2 x' x' have significant difference from the 

pure px 1 ,x 2 ,x[,x' 2 (t) ■ 

We assume the initial state of the external field is in 
a coherence state, \<p) — \a) (noting the condition n < 1 
in the above section, we assume \a\ < 0.5 here). By a 
straightforward calculation, we obtain 

\<j>(t,Xi,X 2 , X[, A' 2 )) 

_ e »*(t,Ai,A 2 ,Ai,Ai) x 



T = 



a + * (A' 1+ A 2 ) 

where the time-varing phase 
"5 



-iuit 



--(A1 + A2) 



(48) 



$ = 



(A' 1+ A' 2 ) 



Ima* 



— (Ai + A2) Im(aV wt ) 



(^) 2 (A' 1+ A' 2 )(A 1+ A 2 ) 



smut 



(49) 



depends on the parameters A'^A'j, Ai and X 2 and 
I [a + 3- (\[ + A' 2 )] e~ iult - 3- (Ai + A 2 )) is also a coher- 
ent state. 

Considering the inner product of coherent state satis- 
fies 

(a \a+p) = exp — \p\ 2 exp [zlm(pa*)] , (50) 
we can express the decoherence factor 

F = exp [iQ] ■ exp [— |T| 2 j (51) 



in terms of the phase factor £1 and its norm exp 
Here, the explicit expression of Q and T can be obtained: 



n = 2 



(Mi +M2 - Mi -M2) 



Ima* 



+2 



— (A'j + A 2 — Ai — A 2 ) Im(a*e 1 ^) 



+ (5) (Mi + M2 - Mi - M 2 ) 
x (Aj + A 2 + Ai + A 2 ) sinwt 



(52) 



and 



-(M1+M2-M1-M2) 
-l(X 1 + X 2 -X' 1 -X' 2 ). 



However, in the case with 



101 << 1,|T| 2 << 1 



(53) 



(54) 



we have (<j) (t; A' l7 A' 2 ; fj,[, fi' 2 ) \<j> (t, Ai, A 2 ; Mi,M2)) - 1 and 
the decoherence effect can be neglected. 

It is observed from the explicit expression for the de- 
coherence factor that O and T are all periodic functions 
of t with frequency u). The amplitudes of ft and T do 
not vary with time (see Fig. 4). Therefore, in case the 
amplitudes of £1 and T are much smaller than 1, the in- 
equalities in Eq. (54) can be satisfied at any instance. 



ReF 




0.95 



FIG. 4. The real part of the decoherence factor Ref as a 
function of t. The unit of t is lj^ 1 . We assume a — and 
H 1 + n 2 - hi - fi 2 ' = Xi + \ 2 - A 1 ' - A 2 ' = 1. The ratio 
^ is denoted as k. It can be seen from the figure that the 
decoherence factor is a periodic function of time. 

Because every Ai or A 2 may be ±1, the sufficient con- 
2 1 1 2 

dition of |T| << 1 is just 8^ << 1 and if ^ is small 

enough so that this condition can be satisfied, then the 

last term of the expression of can be neglected and the 

sufficient condition of 101 << 1 becomes 



16- a « 1. 



(55) 



This is a restriction for the strength of the external field. 
Therefore, we know that if the ratio — is small enough 
(for example, — ~ 0.01) and the external field is weak 
enough (e.g. \a\ « 1), the atomic state decoherence 
caused by the external field can be omitted and the two 
qubit system can evolve under the control of the effective 
Hamiltonian h. 



V. CONCLUSION WITH REMARKS 

In this paper, we have proposed a new protocol for the 
creation of quantum entanglement between two qubits 
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interacting with a high frequency data bus. Unlike the 
previous schemes, our protocal does not require the rota- 
tion wave approximation for the interaction Hamiltonian. 
We can first assume the external field is prepared in the 
coherence state \a). If the frequency of the external field 
is large enough and the strength of it is weak enough so 
that the following conditions 
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uj a ~g«oj, 64^ «1, 
\a\ < 0.5, — \a\ « 1 



(56) 



can be satisfied, an effective coupling between the two 
qubit system can be obtained under the BO approxima- 
tion and the decoherence caused by the cavity field can 
be neglected. 

The explicit expression of the Hamiltonian of the com- 
posite system is given in Eq. (3). If the conditions by Eq. 
(56) are satisfied, the evolution of the two qubit system 
is governed by the effective Hamiltonian 



h = uj a (a 



(57) 



In our previous calculation, we have assumed u> a <~ 
g « uj. This is the sufficient condition to obtain the 
effective Hamiltonian (57). It is easy to prove that if u a is 
smaller than g and g « ui, the effective Hamiltonian 
(57) can be also obtained with the deduction in the above 
sections. To realize a logic gate with high efficiency, we 
suppose our proposal work under the condition u> a << 

2 

g and g « uj so that the interaction strength 2^- is 
comparable with the strength cu a of the free Hamiltonian. 

As we have shown, the decoherence effect is determined 
by the factor fl and |Y| which are periodic functions of 
time. Then if the conditions in Eq. (56) are satisfied, the 
decoherence effect can be neglected no matter how long 
the operation takes. 

The scheme proposed in this paper can be generalized 
to the system that N qubits coupled to a data bus. With 
the same discussions as above, it is apparently that in the 
large detuning and weak coupling case where u) a « uj, 
g « ui, and the number of the qubits is not very large 
so that the condition 2Nuo a « uj is satisfied, an effective 
coupling between the qubits 



can be obtained under the BO approximation. It is 
pointed out that with the "Jj interaction" between N 
qubits, the A-qubit GHZ state can be created easily ( 
[16]). Then a physical qubit of the form a |0) + |1) can 
be encoded into a logical qubit of the form a |000. ..) + 
1 1 1 1 . . .) ( [17]). It is also proved that with the interac- 
tion of this type, Shor code for error correction can be 
easily realized. 



APPENDIX A: THE CALCULATION OF THE 
OVERLAP FACTOR <M(Ai,A 2 ) \N (Ai, A' 2 )) 

In this appendix, we calculate the explicit expression 
of (m(Ai,A 2 ) ^(A^A^)). It can be seen from Eq. (6) 
that the Hamiltonian Hf (Ai, A 2 ) is just the Hamiltonian 
of a displaced harmonic oscillator and can be expressed 
as 



H f (A 1; A 2 ) uA\ iM A XlM - 2^\ 1 \ 2 



(AI) 



where the constant term independent on Ai,A 2 is omitted. 
Here. 



Ax lM = a+ - (Ai +A 2 ) 

UJ 



is a 

tor that satisfies 



displaced 



A 



Ai ,A 2 



boson opera- 
= 1. With Eq. (AI), 



the eigen-state of Hf (Ai, A 2 ) can be written as: 



KAi,A 2 )) 



Vn! 



D 



--(Ai+A 2 ) \n) 



in terms of the displacement operator D [— & (Ai 
and the natural Fock state \n) = -7=^" |0). 



(A2) 

+■ A2)] 
Here, 



|0) is the realistic vacuum state satisfies a |0) = 
and |0(Ai,A 2 )) the displaced vacuum state satisfies 
^Ai,A 2 |0 (Ai, A 2 )) = 0. It is easy to proof that |0 (Ai, A 2 )) 
is actually a coherence state which can be defined as 
|0(Ai,A 2 )) = £[-£(Ai + A 2 )] |0>. 
With 

Eq. (A2), the explicit value of (m(Ai,A 2 ) \n(X' l7 \' 2 )) 
can be obtained easily: 



(m(Ai,A 2 ) |n (A' 1; A 2 )) = (m\D 



.y vSAi,A 2 

— 2 -'\' V 



|n) 



= e 



.19 v Al ,A 2 



.1 3 v Al ,A 2 



w Ai,A 2y / 



i' A 2 \n) 

2l+(n-m) 



1=0 



l\(l + n- m)\ 



sjn{n - 1) ... (m - I + \)\J m (m - 1) .. (to - I + 1). 

This is just Eq. (16) in section III. Substituting Eq. (16) 
into Eq. (14) and Eq. (15), we can get the expression of 
(n) and O^, 2 (n) i.e. Eq. (17) and Eq. (18). 
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